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Data-Based  Control  of  a  Free-Free  Beam  in  the  Presence  of  Uncertainty 

Tannen  S.  VanZwieten,  Gregory  M.  Bower  and  Seth  L.  Lacy 


Abstract — Linear  control  development  is  typically  based 
on  deterministic  models  that  approximate  the  system  under 
consideration.  This  approach  neglects  uncertainty  in  the  system 
response.  System  uncertainty  can  arise  from  a  number  of 
sources  including  disturbances,  noise,  unmodeled  dynamics,  and 
nonlinearity.  This  may  result  in  a  reduction  in  performance 
or  even  instability  in  the  closed  loop  system.  The  goal  of  this 
research  is  to  account  for  measured  uncertainty  in  control 
design.  Our  approach  is  to  tune  a  baseline  controller  using  a 
cost  function  that  balances  performance  and  robustness  given 
measured  system  uncertainty.  The  approach  is  demonstrated 
on  a  free-free  beam,  with  the  goal  of  mitigating  the  flexural 
vibration.  A  lumped  mass  model  is  tuned  to  match  the  ex¬ 
perimentally  measured  Frequency  Response  Function  (FRF)  of 
an  experimental  beam.  This  evaluation  model  and  a  reduced 
order  model  are  used  to  approximate  the  beam  dynamics.  The 
baseline  (LQG)  controller  is  designed  around  the  reduced  order 
model  of  the  beam.  This  controller  is  tuned  according  to  the 
proposed  cost  function  using  the  FRF  and  postulated  variance. 
The  cost  function  includes  closed  loop  performance  and  stability 
robustness  metrics.  The  resulting  baseline  and  tuned  controllers 
are  evaluated  on  lumped  mass  models  consistent  with  the 
measured  data  and  uncertainty. 

I.  INTRODUCTION 

Deterministic  models  are  routinely  developed  as  the  basis 
for  linear  control  development  by  approximating  the  system 
under  consideration.  However,  no  model  is  perfect.  Uncer¬ 
tainty  is  introduced  to  the  system  via  disturbances,  noise, 
unmodeled  dynamics,  and  nonlinearity.  This  can  result  in  a 
reduction  in  performance  or  even  instability  in  the  closed 
loop  system. 

System  uncertainty  can  be  particularly  significant  with 
flexible  structures.  Since  flexible  structures  are  infinite  di¬ 
mensional,  implementable  control  designs  must  be  reduced 
order. 

One  example  of  a  large  flexible  structure  is  the  Deployable 
Optical  Telescope  (DOT)  at  the  Air  Force  Research  Labora¬ 
tory  on  Kirtland  Air  Force  Base  (shown  in  Figure  1).  Data 
based  models  were  created  for  this  lightly  damped,  modally 
dense  structure  [1],  Variance  information  was  captured  to 
quantify  the  uncertainty.  DOT’s  optical  requirements  neces¬ 
sitate  the  use  of  a  very  precise  controller,  but  this  cannot  be 
at  the  expense  of  stability. 
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Fig.  1.  AFRL’s  Deployable  Optical  Telescope 


This  paper  develops  the  idea  of  using  measured  FRF  and 
variance  information  to  tune  a  baseline  controller.  The  goal 
of  the  tuning  algorithm  is  to  find  a  controller  that  can  account 
for  measured  system  uncertainty. 

Our  approach  is  to  tune  a  baseline  controller  using  a  cost 
function  that  attempts  to  balance  closed  loop  performance 
and  robustness  given  measured  system  uncertainty.  In  the 
absence  of  uncertainty,  the  H2  cost  function  gives  optimal 
performance.  However,  there  are  no  guaranteed  stability 
margins  for  a  H2  (or  LQG)  controller  [2],  [3],  and  in  the 
presence  of  uncertainty  the  standard  H2  optimal  controller 
has  been  found  not  to  be  sufficiently  robust  for  practical 
applications  [4].  Therefore,  a  Nyquist-criteria  based  cost 
function  will  be  included  with  the  H2  cost  to  form  the 
total  cost.  A  similar  technique  was  developed  in  [5].  The 
controller  will  be  optimized  based  on  the  FRF  and  postulated 
variance  information.  The  inclusion  of  FRF  variance  allows 
the  stability  and  performance  of  the  closed  loop  system 
to  be  evaluated  in  terms  of  the  probability  of  meeting  a 
performance  goal  and  the  probability  of  being  stable. 

Our  analysis  is  demonstrated  on  a  free-free  beam,  with  the 
goal  of  mitigating  the  flexural  vibration.  Similar  to  DOT,  the 
free-free  beam  contains  an  infinite  number  of  states,  forcing 
the  controller  design  to  be  based  on  a  ROM.  The  controller 
development  will  be  shown  here  for  the  SISO  case,  but  the 
Generalized  Nyquist  Stability  Criterion  [6]  would  allow  it  to 
be  extended  for  MIMO  systems. 

This  paper  begins  in  Section  II  with  a  description  of 
the  experimental  setup,  including  the  spatial  model,  data 
collection  and  validation.  The  Frequency  Response  Function 


collected  from  the  experiment  is  used  to  tune  a  lumped 
mass  model  that  is  developed  in  Section  III.  A  ROM  of  the 
lumped  mass  model  is  also  presented  in  Section  III.  This  is 
used  to  approximate  the  beam  dynamics  and  serves  as  the 
basis  for  the  development  of  the  baseline  (LQG)  controller 
in  Section  IV.  Section  V  introduces  postulated  uncertainty 
that  is  consistent  with  the  unperturbed  and  perturbed  beam 
models.  In  Section  VI  the  controller  design  is  presented  in 
terms  of  cost  functions  and  their  corresponding  sensitivities. 
Finally,  numerical  simulations  are  presented  in  Section  VII, 
followed  by  conclusions  and  future  work  in  Section  VIII. 

II.  EXPERIMENT 

A.  Spatial  Model 

A  spatial  model  describes  the  physical  characteristics 
of  a  structure,  and  serves  as  the  basis  for  the  theoretical 
model  development  [7].  Here,  a  free-free  beam  is  chosen  for 
analysis,  with  a  control  input  applied  at  x  =  a,  as  shown 
in  Figure  3.  The  physical  characteristics  of  this  structure  are 


Fig.  2.  Free-Free  Beam  with  Force  Applied  at  x  =  a 

given  below. 

Dimensions 

length:  l  =  1.27  [m] 
width:  w  =  0.0508  [m] 
thickness:  t  =  0.0127  [m] 

Material  properties  (aluminum) 

Young’s  Modulus:  E  =  6.7569el0  [Pa] 
density:  p  =  2767.99  [kg/m* * 3] 

B.  Setup 

The  beam  was  suspended  from  the  ceiling  like  a  pendulum 
to  simulate  free-free  boundary  conditions,  as  shown  in  Figure 

3.  Control  input  was  applied  via  a  5  lbf  electrodynamic 
shaker  and  a  load  cell  was  used  to  measure  this  input  force. 

The  output  measurement  point  was  collocated  with  the 
driving  point.  The  velocity/force  response  was  measured  with 
a  laser  vibrometer. 

C.  Driving  and  Measuring  Point  Location 

The  location  of  the  driving/measuring  point  was  strate¬ 
gically  chosen  based  on  the  first  four  mode  shapes.  These 
shapes,  calculated  from  Euler-Bernoulli  free  vibration  solu¬ 
tion  [8],  are  shown  in  Figure  4  and  Figure  5.  The  small 
green  circle  denotes  the  chosen  I/O  location,  a  =  21.27  cm. 
This  point  was  chosen  near  the  peak  of  third  and  fourth 
modes  so  that  the  magnitude  of  vibration  for  these  modes 
was  comparable  to  that  of  the  first  and  second  modes  (which 
is  inherently  larger).  This  choice  facilitates  the  excitation  of 


Fig.  3.  Experimental  Setup 


Mode  Shapes 


Fig.  4.  Mode  Shapes:  First  and  Second  Modes 


the  smaller  modes  in  such  a  way  that  the  output  response 
will  not  be  dominated  by  larger  amplitude  modes. 

D.  Frequency  Response  Function 

Three  sets  of  frequency  response  functions  were  collected 
with  the  following  frequency  intervals  and  bandwidths, 


A/i 

=  0.25  Hz 

BWi  :  0  -  800  Hz, 

A/2 

=  0.125  Hz 

BW2  :  0  -  400  Hz, 

A/a 

=  0.0625  Hz 

BW3  :  0  -  200  Hz. 

The  mobility  (velocity/force  FRF)  collected  experimen¬ 
tally  was  averaged  over  all  three  data  sets,  as  shown  in 
Figure  6.  This  was  integrated  to  give  the  compliance  (dis¬ 
placement/force).  All  future  calculations  are  based  on  the 
average  FRF  values  for  compliance. 


Mode  Shapes 


distance  along  beam,  x  [m] 


distance  along  beam,  x  [m] 


Fig.  5.  Mode  Shapes:  Third  and  Fourth  Modes 


Theoretical  natural  frequencies  [8]  matched  experiment  val¬ 
ues  for  the  first  four  modes  with  less  than  one  percent  error. 

III.  LUMPED  MASS  MODEL 

The  lumped  mass  model  development  relies  on  the  spatial 
model  of  the  beam  given  in  Section  II- A.  This  model  gives 
transient  response  characteristics  of  the  beam  in  addition  to 
the  modal  model  (steady  state). 

For  a  lumped  mass  model,  the  beam  is  split  into  ne 
elements  with  nn  =  ne  +  1  nodes,  as  shown  in  Figure  7. 
Each  node  is  defined  by  one  translational  and  one  rotational 
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Fig.  6.  FRF  for  Free-Free  Beam 


Fig.  7.  Elements  and  Nodes  for  a  Lumped  Mass  Model 


degree  of  freedom,  giving  the  entire  beam  2 nn  degrees  of 
freedom  and  Tnr,  states.  Each  element  is  defined  to  have  a 
length  le  = 


A.  Mass  Matrix 


A  simple  lumped  mass  matrix  is  chosen,  where  the  mass 
of  each  element  is  chosen  in  proportion  to  the  number  of 
elements  in  the  model.  This  is  known  as  an  inconsistent- 
mass  matrix,  since  it  is  not  fashioned  from  a  set  of  shaping 
and  displacement  functions  such  that  it  is  consistent  with  the 
stiffness  matrix  [9],  The  mass  applied  to  the  node  on  the  end 
is  given  by 


?7)end  _  P-Ale 
lne  ~  r, 


(2) 


corresponding  to  half  the  mass  of  the  adjoining  element.  For 
all  other  elements  the  mass  is  given  by 


E.  Validation 

For  any  beam  there  will  be  an  infinite  number  of  nor¬ 
mal  modes,  where  each  mode  is  associated  with  a  natural 
frequency.  The  mode  of  lowest  frequency  is  found  from 
the  experimental  compliance  to  occur  at  1.5  Hz  and  is  a 
rigid  body  mode  rather  than  a  flexural  mode.  This  mode 
is  associated  with  the  support  conditions,  and  in  this  case 
corresponds  to  a  pendulum-like  motion  of  the  beam.  For  our 
analysis,  the  first  four  dominant  flexural  modes  are  selected 
from  the  FRF. 

/i  =  39.75  Hz, 

f-2  =  109.39  Hz, 

H  =  213.5  Hz,  (1) 

fA  =  352.0  Hz. 

These  frequencies  correspond  to  points  of  maximimum 
response  of  the  average  FRF.  The  sharp  peaks  at  these 
frequencies  indicate  that  the  beam  is  very  lightly  damped. 


me  =  pAle, 


(3) 


since  there  are  two  adjoining  elements.  The  mass  moment  of 
inertia  for  each  element  about  the  node  is  found  by  taking  the 
inertia  of  half  of  the  adjoining  element(s)  about  the  node.  In 
the  case  of  the  end  nodes,  where  the  node  connects  to  only 
one  element,  we  find 


rend  _  pAle 

e  ~  24 


(4) 


For  all  other  nodes,  the  inertia  corresponds  to  that  of  two 
half  elements  taken  about  the  node  of  interest,  and  hence 
the  above  quantity  is  doubled, 


Ie  = 


pAll 

12 


(5) 


The  mass/inertia  matrix  for  the  entire  beam  is  given  by 


M  =  diag  [mfd  7eend  me  Ie  ...  me  Ie  mfd  7eend]  ,  (6) 
where  M  e  R2"=x2n=. 


B.  Stiffness  Matrix 

The  stiffness  assigned  to  each  nodal  DOF  corresponds 
to  the  stiffness  of  the  element(s)  neighboring  the  node  of 
interest.  Therefore  the  stiffness  of  each  element  must  be 
characterized  and  then  related  to  the  nodal  linear  and  angular 
displacement. 

The  forces  and  moments  are  applied  at  the  nodes  on  either 
end  of  the  element,  xe  =  0  and  xe  =  le,  as  shown  in  Figure 
8.  The  forces  and  moments  at  nodes  1  and  2  are  related  to 


Fig.  8.  Forces  and  Moments  on  a  Beam  Element 


as  shown  in  [11]. 

For  a  beam  that  contains  multiple  elements,  the  stiffnesses 
for  each  element  given  by  (8)  must  be  assembled  such 
that  the  nodes  of  the  entire  beam  are  represented  in  a 
single  matrix.  Note  that  the  first  two  rows/columns  of  I\e 
correspond  to  the  linear  and  angular  stiffness  of  node  1,  while 
the  last  two  rows/columns  correspond  to  node  2.  Also  note 
that  node  2  of  an  element  represents  the  same  node  as  node 
1  of  the  adjoining  element,  provided  that  there  is  an  element 
attached  to  either  end  (i.e„  the  node  considered  does  not 
represent  a  boundary  condition). 

The  block  diagonal  matrices  of  the  stiffness  for  each  node 
with  zeros  at  all  other  nodes  are  summed  over  all  nodes  to 
give  the  stiffness  matrix  for  the  entire  beam. 

C.  Equation  of  Motion 

The  equation  of  motion  for  the  lumped-mass  beam  is 


the  stiffness  of  the  element  by 


h  ' 

12  6  le  -12  6  le 

di 

mi 

El 

6  le  All  — 6/e  2/2 

01 

h 

m2 

-12  -6  le  12  -64 

_  64  2 /2  -64  Al\ 

<4 

02 

(7) 


The  shape  functions  corresponding  to  (7)  illustrate  the 
effect  of  nodal  displacements  and  rotations  on  the  transverse 
displacement  of  a  given  element  along  its  length  [10].  The 
plot  of  the  shape  functions  along  the  length  of  the  beam  is 
shown  in  Figure  9. 


Md(t)  +  Cd(t)  +  Kd(t)  =  F(t ),  (9) 

where  d  =  [  d\  0i  d2  f2  ■■■  dUm  ]  €  R2nm 

is  the  linear  and  angular  displacement  at  each  node  and  F  e 
R2n™  is  the  vector  of  forces  and  moments  applied  at  each 
node.  Defining  x\  =  d,  x2  =  d,  and  x  =  [x\  x2]T ,  this  may 
be  re-written  in  state  space  form  as 

x(t)  =  Ax(t)  +  Bu(t) ,  (10) 

where 


Fig.  9.  Shape  Functions  for  a  Beam  Element 


By  definition,  the  stiffness  is  the  force  divided  by  the 
displacement,  so  from  (7)  we  see  that  the  element  stiffness 
is  simply 
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The  lumped  mass  FRF  converges  to  the  experimental  FRF 
as  the  number  of  nodes  is  increased.  The  model  was  tuned 
by  adjusting  the  damping  coefficients  along  the  modes.  The 
addition  of  damping  results  in  a  stable  open  loop  system, 
which  is  important  for  the  existence  of  the  H2  norm  that 
will  be  used  as  the  performance  metric  in  the  controller 
cost  function.  This  stable  open  loop  assumption  is  also 
an  important  consideration  for  the  Nyquist-based  stability 
metric. 

The  FRF  for  the  lumped  mass  model  with  55  elements 
and  linear  and  angular  damping  coefficients  of  c  =  0.1808 
and  c  =  1.33e  —  6,  resp.  is  shown  in  Figure  11. 

Since  the  focus  for  the  model  development  and  control 
in  this  paper  is  on  the  flexural  modes  of  vibration,  the  rigid 
body  mode  was  eliminated  via  a  high  pass  Butterworth  filter. 
This  is  shown  in  Figure  10  for  the  lumped  mass  model.  The 
same  filter  was  applied  to  the  experimental  data. 


TRANSVERSE  VIBRATION  MODEL 


Fig.  10.  Eliminate  Rigid  Body  Mode  With  a  High  Pass  Filter 


E.  Reduced  Order  Model 

A  ROM  is  developed  to  approximate  the  dynamics  given 
by  the  lumped  mass  model  in  (10).  The  purpose  is  to  reduce 
the  model’s  computational  complexity  as  much  as  possible 
while  retaining  the  major  system  characteristics.  The  baseline 
and  robust  controllers  will  be  developed  around  the  ROM. 

The  lumped  mass  model  was  transformed  into  modal 
form.  The  states  with  a  small  contribution  to  the  input- 
output  response  were  identified,  isolated,  and  eliminated.  The 
FRF  for  a  ROM  of  order  nr  =  12  and  a  lumped  mass 
model  with  55  elements  (220  states)  is  plotted  alongside 
the  beam’s  experimental  FRF  in  Figure  11.  The  SISO  ROM 
is  represented  by  (Ar,  Br,Cr,  Dr)  with  Ar  £  M"rX”r, 
Br  £  R”’’,  Cr  £  R”r,  and  Dr  £  R. 


ROM  FRF  at  x=a 


Fig.  11.  Lumped  Mass  Model.  Reduced  Order  Model  and  Experimental 
Compliances 

IV.  BASELINE  CONTROLLER 

A.  Linear  Quadratic  Gaussian 

An  LQG  controller  may  be  written  in  state  space  form  as 

xc(t)  =  Acxc{t)  +  Bcuc(t)  (12) 

u  =  Ccxc(t)  (13) 


where  A, 

c  e  Rn-xn- 

,  Bc  £  Rn-xl,  Cc  £  R1™-,  Dc 

£  R, 

and 

Ac 

=  Ar-  BrK  -  LCr 

(14) 

Bc 

=  L 

(15) 

Cc 

=  —K. 

(16) 

where  L 

is  the  Kalman  filter  gain  and  K  is  the  LQR 

static 

state  feedback  gain. 

B.  Parameterization 

The  baseline  controller  is  developed  using  the  SISO  ROM 
that  contains  nr  states.  Hence,  the  standard  LQG  form  results 
in  n2  +  nr  parameters  that  will  need  to  be  tuned  using  the 
algorithm  in  the  following  section.  Conversely,  for  the  same 
controller  in  modal  form,  there  exist  3 nr  parameters  to  esti¬ 
mate.  Therefore  by  transforming  the  baseline  controller  into 
modal  form  there  is  a  reduction  of  n2  —  2 nr  parameters.  For 
a  ROM  with  nr  =  12  states,  the  LQG  controller  in  standard 
form  results  in  156  parameters,  but  once  transformed  into 
modal  form  there  are  only  36  parameters  -  a  reduction  of 
120  parameters. 

V.  UNCERTAINTY 

Uncertainty  in  systems  is  caused  by  a  number  of  different 
sources,  including  disturbances,  noise,  unmodeled  dynamics 
and  nonlinearities.  The  effect  of  this  uncertainty  may  be  a 
reduction  in  performance  or  even  instability  of  the  closed 
loop  system. 

We  introduce  uncertainty  in  the  form  of  postulated  vari¬ 
ance  information.  The  variance  is  the  square  of  this  standard 
deviation,  a2.  This  standard  deviation  may  be  used  to  de¬ 
fine  confidence  intervals  around  the  experimental  FRF.  For 
example,  a  95%  confidence  interval  is  given  by  1.96ct  and  a 
99%  confidence  interval  is  given  by  2.58er. 

For  our  purposes,  the  postulated  variance  is  simply  cal¬ 
culated  from  the  compliances  of  unperturbed  and  perturbed 
lumped  mass  models  that  will  be  used  to  test  the  controller 
robustness  in  Section  VII.  In  particular,  1.96cr  was  defined 
such  that  it  smoothly  encompasses  the  deviation  of  the 
perturbed  system,  as  shown  in  Figure  12.  This  corresponds 
to  a  95%  confidence  interval  around  the  displacement  of  the 
unperturbed  system,  shown  in  Figure  13. 

VI.  ROBUST  CONTROLLER  DESIGN 

The  robust  controller  design  considered  here  augments  an 
H-2  performance  cost  with  a  Nyquist-based  stability  cost. 
Both  cost  functions  are  data-based.  hence  the  tuning  relies 
completely  on  the  FRF.  Note  that  for  the  remainder  of  the 
paper  the  specific  dependence  on  juj  or  /:.  where  obvious, 
may  be  omitted  for  simplicity  of  exposition. 

A.  Closed  Loop  System 

Let  the  state  space  representation  of  the  ROM  be  repre¬ 
sented  by, 

xr  =  Axr  +  Buu, 

z  —  Czxr  D zuu ,  (17) 

y  =  CyXr  +  DyuU, 


3.5 


Fig.  12.  Postulated  Standard  Deviation 


frequency,  f  [Hz} 


X 

A 

BuCc 

X 

Xc 

BcCy 

(Ac+BcDyuCc) 

Xc 

z  =  [Cz  DzuCc  ] 
y  =  [Cy  DyuCc  ] 


X 

ir¬ 


respectively.  The  above  closed  loop  matrices  are  denoted  by 

A^cl\  d  ,  and  C^cl). 

B.  Cost  Functions 

The  data-based  performance  cost  is  taken  as  the  weighted 
sum  of  the  //2  for  the  output  and  performance  variables, 


where 


Jp  OLyJpy  T  azJpz, 


len(u. 


Jpv=~  55  tr  {HyU^n^U^k)}  Aw, 

k= 1 
len(co) 


jPx  =  ^  55  tr {Hx(ju)Hx(juk)}  Aw, 


(23) 

(24) 

(25) 


fc=i 


and 

ByijiOkj  —  Gyu(ju>k)  S(j Uy'k)  . 


Hz(jujk )  =  Gzu(juk)S(jujk)  (26) 


-  (r  -  Gyu(juk)K(juk))  ,  (27) 

Aw  is  the  frequency  spacing  between  consecutive  points,  and 
H  denotes  the  conjugate  transpose.  The  relative  weighting 
between  az  and  ay  determines  the  emphasis  placed  on  the 
output  verses  the  entire  state  vector.  This  cost  maximizes 
performance,  but  has  no  guaranteed  stability  margin. 

The  Nyquist-based  stability  cost  is  given  by 


Js  =  ~  V  W(wk) 

IT  J 


len(uj) 


l 


fc= i 


d2(ju>k) 


Aw 


where 


d(jvk)  = 


|1  +  GyU(juik)K{juk)  | 


(28) 


(29) 


Fig.  13.  Perturbed  and  Unperturbed  Systems  with  95%  Confidence  Interval 


where  xr  denotes  the  state  of  the  ROM  and  z  is  a  per¬ 
formance  measure.  The  corresponding  transfer  functions  are 
given  by 

GyU(s)  =  Cy  (si  —  A)-1  Bu  +  Dyu,  (18) 

Gzu(s)  =  Cz{sl~  A)-1  Bu  +  Dzu.  (19) 

Combing  the  above  ROM  plant  equations  from  (17)  with 
the  controller  given  by  (12)  and  (13),  the  closed  loop  state 
equation  becomes 


(20) 


and  the  closed  loop  performance  and  output  equations  are 


(21) 


r(wfc) 

and  T’(wfe)  corresponds  to  a  95%  confidence  interval  at  the 
frequency  of  interest, 

T(wfc)  =  1.96er(wfc).  (30) 

This  is  a  slightly  modified  version  of  the  Nyquist-based 
cost  function  found  in  [5],  It  takes  into  account  the  stability 
criterion  as  well  as  the  ability  to  handle  model  uncertainty. 

Nyquist  analysis  for  SISO  systems  is  a  graphical  stability 
test  that  is  based  on  a  plot  of  the  loop  gain  G(jw)/\(jw) 
along  a  frequency  contour.  For  closed  loop  stability,  the 
number  of  encirclements  of  the  (—1,0)  critical  point  must 
be  equal  to  the  number  of  open  loop  unstable  poles.  Since 
the  system  under  consideration  is  open  loop  stable,  we 
require  there  to  be  zero  counterclockwise  encirclements  of 
the  critical  point  for  closed  loop  stability.  The  distance  metric 
(29)  is  simply  the  magnitude  of  the  distance  to  the  critical 
point  divided  by  the  confidence  interval.  Taking  a  weighted 
sum  of  the  inverse  of  this  distance  metric  squared  means  that 
the  cost  function  will  decrease  as  the  curve  moves  farther 
from  the  critical  point  and/or  has  a  smaller  variance. 

The  total  cost  is  a  weighted  sum  of  the  stability  and 
performance  costs, 

J  —  OiJp  PJS. 


X, 


(22) 


(31) 


Equation  (31)  quantifies  the  stability  verses  robustness  con¬ 
trol  problem.  The  relative  weighting  between  a  and  (3  deter¬ 
mines  the  emphasis  placed  on  performance  verses  stability 
robustness.  The  beam’s  controller  will  be  tuned  with  this  cost 
function  using  the  beam’s  FRF  and  postulated  variance  data. 


The  Nyquist  stability  plot  (see  Figure  14).  however,  reveals 
that  the  system  becomes  unstable  when  a  multiplicative 
perturbation  of  0.93  is  applied  to  the  stiffness  term.  This 
is  the  same  perturbation  as  was  used  to  characterize  the 
uncertainty  in  Section  V. 


C.  Sensitivity 

The  sensitivity  is  determined  by  taking  the  gradient  of 
the  cost  functions  with  respect  to  the  controller  parameters. 
The  output  performance  sensitivity  for  an  arbitrary  controller 
variable,  pi,  is  given  by 


f)  I 

VUp'i 

dpt 


len{uj) 


=  -  E  tr  Re 


k= 1 


Hy  (jujk)  }  Aw, 


(32) 


where 


dJk^G  m) 

dPi  GyudPl ’  (33) 

=  SGyu^S,  (34) 

dpi  dpi 

0K  =  ^$cBc  +  Cc<f>c^<I>cBc  +  Cc<I>c^,(35) 

OPi  OPi  OPi  OPi 

$c  =  (si -Ac)-1.  (36) 


dJ, 


is  similarly  defined.  The  sensitivity  of  the  stability  cost 
with  respect  to  pi  is 


len(u) 

aA±  =  iy  w1 

OPi  7 T  t—1 


1 


k= 1 


d 2  |1  +  GyuI\\ 


OK 

OPi 


Gyuyyy  Aw,  (37) 


where  js  defined  in  (35).  This  is  a  modified  version  of 
the  sensitivity  given  by  [5]. 

The  complete  sensitivity  is  given  by 


dJ_  _ 

dpi 


+  ay 


8Jpy\ 

dpi  ) 


(38) 


Each  controller  parameter,  p,.  is  modified  according  to  the 
above  total  sensitivity  functions  such  that  the  total  cost 
function  is  locally  minimized. 


VII.  NUMERICAL  SIMULATIONS 

Numerical  simulations  were  performed  on  a  lumped  mass 
model  with  55  elements.  The  input  and  output  were  non- 
colocated;  the  input  was  at  node  10  (corresponding  to  the 
experimental  location),  and  the  output  displacement  was  cho¬ 
sen  to  be  at  node  26.  The  baseline  controller  was  developed 
around  a  ROM  with  12  states  and  the  parameters  from  this 
controller  provided  the  initial  values  for  the  controller  tuning. 
Nyquist  plots  for  both  controllers  are  shown  for  the  nominal 
and  perturbed  plants. 

A.  Baseline  Controller 

The  optimization  gains  for  the  LQR  controller  were  set 
to  be  Q  =  5/i2xi2  and  R  =  le  —  12.  The  Kalman  filter 
process  and  measurement  noise  covariances  were  set  to 
P  =  10  and  M  =  1,  respectively.  The  closed  loop  response 
for  the  nominal  system  results  in  excellent  performance. 


Nyquist  Diagram 


Real  Axis 


Fig.  14.  Nyquist  Plot  for  Baseline  Controller  with  Perturbed  and  Unper¬ 
turbed  System 


B.  Robust  Controller 

The  performance  measure,  z,  was  taken  to  be  the  state  x, 
and  therefore  Cz  =  I  and  Dzu  =  0.  The  weighting  W(wk) 
was  taken  to  be  unity  when  the  real  component  of  the  loop 
gain  was  positive,  and  20  when  it  was  less  than  or  equal  to 
zero.  Furthermore,  the  individual  cost  function  weights  were 
set  to  az  =  1,  av  =  lelO,  a  =  1.  and  (3  =  le4. 

From  the  Nyquist  plot  (Figure  15)  we  can  see  that  the 
tuned  controller  regains  closed  loop  stability  in  the  presence 
of  uncertainty. 


Nyquist  Diagram 


Fig.  15.  Nyquist  Plot  for  Tuned  Controller  with  Perturbed  and  Unperturbed 
System 
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J- total 

Jp-  performance 

•/^-stability 

Baseline 

Tuned 

3.3199e6 
2.94 13e6 

331.9889 

294.1340 

8.9804e-8 

3.8619e-8 

The  tuned  controller  gave  a  reduction  in  both  the  perfor¬ 
mance  and  stability  costs. 


VIII.  CONCLUSION  AND  FUTURE  WORK 

In  conclusion,  a  baseline  controller  was  tuned  using  a  data- 
based  cost  function  that  includes  a  Nyquist-based  stability 
robustness  criterion  with  the  more  standard  Ifj  closed  loop 
performance  cost.  This  used  FRF  and  variance  data  for  the 
nominal  plant  to  tune  the  controller.  The  tuned  controller 
was  applied  to  both  perturbed  and  unperturbed  lumped  mass 
models  of  a  free-free  beam.  The  results  showed  that  the 
closed  loop  system  using  the  baseline  controller  became 
unstable  in  the  presence  of  uncertain  dynamics.  However, 
the  addition  of  the  Nyquist-based  cost  for  the  controller 
optimization  forced  the  closed  loop  system  to  remain  stable 
for  perturbations  within  the  95%  confidence  interval. 

For  future  work,  the  authors  plan  to  determine  the  response 
characteristics  for  various  systems  within  the  family  of 
allowable  plants  (corresponding  to  the  uncertainty).  This  will 
include  various  initial  conditions  as  well  as  input  and  output 
disturbances.  Furthermore,  this  work  could  be  expanded  and 
implemented  for  MIMO  systems.  It  would  also  be  interesting 
to  apply  the  control  theory  to  more  complex  flexible  struc¬ 
tures,  such  as  the  Deployable  Optical  Telescope. 
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C.  Comparison 

The  total,  performance  and  robustness  cost  functions  are 
shown  in  Table  1  for  both  the  baseline  and  tuned  controllers. 

Table  1:  Cost  for  Baseline  and  Tuned  Controllers 


